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Abstract 

In the experiment eight populations of Picea abies were chosen at 100 m intervals between 500 m and 1200 m alti- 
tude a.s.l.. In each population wood core samples were collected from 14-19 trees (126 cores total), and measured 
using a Corim Maxi device. At four of the eight sites (every 200 m in elevation between 500 m and 1100 m a.s.l.), the 
diversity of ground vegetation was evaluated, and temperature was recorded at every 100 m of altitude. 

The highest average radial increment of spruce occurred between the altitudes 800-1000 m a.s.l., which is prob- 
ably the optimum for spruce. The larger increment indices observed at higher altitudes may signify a high growth 
potential of spruce. It may also suggest a recent upward shift of the optimum growth zone for this tree species. 

In 15 phytosociological records, the presence of 148 plant species forming plant associations: Dentario glan- 
dulosae-Fagetum typicum (sub-mountainous and mountainous form) and Abieti-Piceetum, and community Abies 
alba-Rubus hirtus, was documented. No relationship was found between ground vegetation species diversity 
(expressed by Shannon-Wiener index) and levels of stand diversity. The vegetation species diversity varied with 
the elevation above sea level: the highest plant diversity was found at 500 m a.s.l., and decreased with increasing 
altitude. 

The potential increase in air temperatures may result in changes to the altitudinal range of many plant species 
including trees, and consequently in an upward shift of the boundaries of plant zones; in this case the sub-mountain- 
ous and lower mountainous forest zone. In this region, the optimal zone for Norway spruce may be restricted to the 
highest elevations. 
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Introduction 

The climate change forecasted could lead to a transfor- 
mation of forest vegetation and related organisms, e.g. 
herbivorous insects and their natural enemies. Moun- 



tain vegetation appears to be particularly vulnerable to 
such changes (Hattenschwiler and Korner 1995; Penue- 
las and Boada 2003). The main purpose of the BAC- 
CARA project, conducted within the European Com- 
mision 7. Framework Programme, was to develop the 
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tools that will enable forest managers and policy mak- 
ers to evaluate the risk of European forest biodiversity 
and productivity loss under climate change. One of the 
experiments was established in Poland (Carpathian 
Mountains, Radziejowa Massif), where the influence 
of variable climatic conditions on forest growth, stabil- 
ity and resistance, was assessed (Ambrozy and Grodzki 
2013). This paper addresses the results concerning only 
two aspects of the study: vegetation diversity, and the 
radial increment of trees. 

The growth of trees is determined primarily by site 
conditions, but usually influenced also by their health 
status and vitality. A decrease in annual radial incre- 
ment is one of the first symptoms of weakening in trees, 
before others such as crown discoloration become vis- 
ible. For this reason the increment tendency (calculated 
as a quotient or index) can be used to assess tree vitality 
and/or its reaction to stressors, under the condition that 
data is collected from dominant trees, in order to avoid 
the effects of competition. Such a method has previously 
been used to assess the effects of air pollution (Zawada 
1983), or climate variation (Zawada 1994) on trees. 

Between 1979-1984 the Norway spruce stands in 
the upper forest zone on Radziejowa were damaged by 
insects of the genus Cephalcia Panzer (Hymenoptera: 
Pamphiliidae), mainly C. alpina Klug (Jachym 2007). 
In some stands that survived the outbreak with me- 
dium or light defoliation, an assessment of increment 
response was still possible. 

In this paper the authors aim to describe the plant 
diversity and diameter at breast height (DBH) increment 
of spruce trees, in relation to changes in climatic condi- 
tions that correspond to altitude, in order to contribute 
to the continued evaluation of potential forest responses 
to predicted climate change. 

Material and methods 

The field studies were carried out in the years 2010-2012. 
The study area was located in the Carpathians in the 
Radziejowa Massif (Beskid Sqdecki, Poland). 

The Radziejowa Massif was formed from the rocks 
of the Carpathian flysch. It is a sequence of shales regu- 
larly interbedded with sandstones (Margielewski 2000). 
Rocks of the Carpathian flysch make up the bedrocks of 
clay soils. According to the climatic conditions, brown 



soils develop at low and medium elevations, and pod- 
zolic soils at the highest elevations (Maciaszek 2000). 
The annual precipitation in the Beskid Sqdecki region 
is approximately 800-1100 mm, and is a function of 
elevation (being lowest at 500 m a.s.l., and highest at 
1200 m a.s.l.). Mean annual temperature is approxi- 
mately 4-7°C, and is negatively correlated with eleva- 
tion (Brzezniak and Czemerda 2000). 

The study area was located in two vegetation 
stages in different altitudinal zones. The lower forest 
stage (dominated by beech) extends up to an altitude of 
1100 m a.s.l., while the upper forest stage (dominated by 
Norway spruce) on the southern slopes of Radziejowa 
Massif, is well developed above 1100 m a.s.l. (Staszk- 
iewicz 2000). 

The study area is covered by mixed Norway spruce 
Picea abies (L.) H. Karst., European beech Fagus syl- 
vatica L., and silver fir Abies alba Mill, stands, but 
Norway spruce is the dominant species. The main for- 
est management strategy has combined objectives: 
predominantly to protect soil and water quality, and 
to promote ‘close-to-nature forestry’ (Duncker et al. 
2012). These strategies require maintenance of the for- 
est cover (i.e. ‘continuous cover forestry’). This means 
that in most cases stands are naturally regenerated un- 
der a canopy of older trees (Szewczyk 2000). 

Eight populations of Picea abies growing on el- 
evations from 500 m to 1200 m a.s.l. (at every 100 m of 
elevation) were sampled and measured (fig. 1, tab. 1). 
In order to obtain the largest time-series possible, the 
oldest spruce stands were chosen from the limited 
number available at each elevation. In each population 
a single wood core was collected from each tree us- 
ing a Pressler borer at 1.3 m above ground level, from 
14-19 trees with an above-average DBH, in order to 
determine the response of spruce growth to various 
climate conditions. The cores were measured using 
a Corim Maxi device. At each altitude the average 
30-year increment of Norway spruce was calculated 
along with increment indices, based on the compari- 
son of mean yearly increment within selected decades 
(Zawada 1983, modified). Two indices were calculated 
at all altitudes: WP00 (a ratio between 1991-2000 and 
1981-1990), and WP\0 (a ratio between 2001-2010 and 
1991-2000), in order to determine the increment ten- 
dency over last two decades. The indices were calcu- 
lated as follows: 
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Fig. 1 . Location of core sampling in Radziejowa Massif 
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WPOO = 



WP d 



^1976-1985 

^1966-1975 



where: 

/| 99 , 2000 ~ diameter increment in the decade 1991-2000 
[mm], 

^1981 1990 — diameter increment in the decade 1981-1990 
[mm]. 

JJ/piQ- ^2001-2010 

A 99 1-2000 

where: 

/ 2001-2010 ~ diameter increment in the decade 2001-2010 
[mm], 

-^i 99 i_ 2 ooo — diameter increment in the decade 1991-2000 
[mm]. 

Additionally, for the samples taken in the upper for- 
est stage (1100 and 1200 m a.s.l.), indices were calcu- 
lated to correspond with the timing of Cephalcia sp.: 
“damage” - WP d (1976-1985 vs. 1966-1975), and “re- 
covery” - WP r (1986-1995 vs. 1976-1985), in order to 
assess the reduction in increment due to defoliation, and 
the recovery process over the following years. These in- 
dices were calculated as follows: 



where: 

A 976 - 1985 - diameter increment in the decade 1976-1985 
[mm], 

/ ] 9 f ,6 , 975 - diameter increment in the decade 1966-1975 
[mm]. 

yyp _ ^1986-1995 
r / 

1976-1985 

where: 

A 986 - 1995 - diameter increment in the decade 1986-1995 
[mm], 

A 976 - 1985 — diameter increment in the decade 1976-1985 
[mm]. 

Level of defoliation was determined using archive 
data (Jachym 2008). 

Temperature loggers were installed at each eleva- 
tion, in order to determine the annual average tempera- 
ture variation as it relates to the altitude a.s.l. (Urban 
2010 ). 

At four altitudes between 500 m to 1100 m a.s.l. 
(every 200 m), the diversity of ground vegetation 



Tab. 1 . Characteristics of sampled spruce stands 



Com- 

partment 


Coordinates 


Altitude 

a.s.l. 

[m] 


Prop, 
of spruce 
[%] 


Mean age 

of stand from FMP* 
[years] 


Stand 

density 


Number of 
sampled trees 


Mean DBH 
of sampled 
trees [cm] 


Mean age 
of sampled trees 
[years] 


97 a 


49°25’27”N 

20°28’37”E 


500 


40 


60 


0.6 


15 


46.8 


98 


36 f 


49°24’54”N 

20°33’08”E 


600 


30 


45 


0.7 


14 


46.0 


61 


36 a 


49°25’25”N 

20°33’32”E 


700 


40 


regeneration phase 


0.4 


15 


41.6 


83 


44 a 


49°26’45”N 

20°33’32”E 


800 


30 


50 


0.7 


15 


41.8 


69 


28 b 


49°26’22”N 

20°34’10”E 


900 


60 


65 


0.7 


19 


41.9 


62 


19 d 


49°26’53”N 

20°35’11”E 


1000 


20 


45 


0.7 


14 


39.9 


81 


17a 


49°27’05”N 

20°35’29”E 


1100 


60 


regeneration phase 


0.2 


17 


45.7 


108 


14 d 


49°27’03”N 

20°35’56”E 


1200 


80 


60 


0.9 


17 


32.6 


65 



4 FMP - Forest Management Plan. 
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was evaluated using a set of 15 experimental plots 
(20 x 20 m squares, 0.04 ha). At each elevation four 
plots were selected according to a tree species composi- 
tion gradient: from pure/nearly pure Norway spruce to 
mixed (four diversity levels per altitude, with the excep- 
tion of 1100 m a.s.l., which had three due to limited tree 
species diversity). This means, that at each elevation the 
first plot represented pure (or nearly pure) spruce stand 
and the last plot - stands with the richest composition 
(up to seven species). 

On each plot the vegetation was described via phy- 
tosociological record, using the Braun-Blanquet meth- 
od, in order to determine plant associations and com- 
munities. For each phytosociological record, the Shan- 
non-Wiener (//’) species diversity index was calculated 
using the formula as follows: 

H' = -Y,p i lnp i 

where: 

Pi is the proportion of individuals of species i. 

The species coverage values in the five degree 
Braun-Blanquet scale were converted into the van der 
Maarel (1979) scale. 

Data was analyzed via correlations, one-way ANO- 
VA and post-hoc tests (Tukey RIR) using STATISTICA 
5 software (StatSoft 1997). 

Results 

Mean annual temperature, estimated from the data col- 
lected in the period 15.05.2010 to 14.05.2011 for the 
Radziejowa Massif, is negatively correlated with alti- 
tude (r 2 = 0.94, p < 0.001) (fig. 2). The highest mean 
value (6.92°C) was found at 500 m a.s.l., decreasing 
constantly with increasing altitude to the lowest mean 
value (3.71°C) at 1200 m a.s.l. On the points located in 
the upper forest stage (1100-1200 m a.s.l.), the mean 
temperature (4.04°C) was much lower than in the zones 
500-700 m (6.52°C), and 800-1000 m (5.37°C) m a.s.l., 
and the differences between zones were statistically 
significant (ANOVA F = 34650.83; p < 0.001). 

In total 126 spruce trees were sampled. The DBF! of 
sampled trees was 41.9 cm on average, varying between 
individual locations within a range 32.6-46.8 cm. The 
mean age of these trees varied from 61 to 108 years old, 



and according to the Forest Management Plan (tab. 1), 
was generally higher than the stand age. 




Fig. 2. Thermal patterns - mean annual temperature at the 
elevations 500-1200 m a.s.l. 




Fig. 3. Average 30-year radial increment of Norway spruce 
at individual altitudes 



The 30-year average DBH increment of eight spruce 
populations studied significantly differed between indi- 
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vidual altitudes (ANOVA F = 3.586; p < 0.01), being 
highest at the altitudes 800 and 1000 m above sea level 
(only these two altitudes were statistically different), 
while the stands at other elevations had a much lower 
average increment. The trend line of the radial incre- 
ment is at its maximum between the elevations where 
stands were characterized by highest growth (fig. 3). 

The increment index calculated using the data from 
last three decades (WP00 and JVP10) has the highest 
value in the stands growing at the highest locations, at 
more than 900 meters above sea level (tab. 2). These 
stands are characterized by increasing growth in recent 
decades, while the others have decreased in growth or 
displayed significant variation in the subsequent dec- 
ades. Generally, the yearly increment did not differ 
significantly between decades at the locality level, re- 
gardless of the value of calculated indices, except in two 
cases where differences at p < 0.05 were found (tab. 2). 

Damage caused by Cephalcia during the outbreak 
from 1979-1984 affected only the stands in the higher 
forest zone. Among those sampled, only the stand grow- 
ing at 1200 m a.s.l. was affected by moderate (31-60%) 
defoliation, resulting in a reduced radial increment be- 
tween 1979-1981, which was especially visible when 
compared with an undamaged stand (fig. 4). The differ- 



ences between the decades analyzed were statistically 
significant only in the damaged stand, which demon- 
strated rapid recovery in increment during the decade 
following defoliation (tab. 3). 

Tab. 2. Values of radial increment indices for the stands on 
the observation plots in the Radziejowa, and the significance 
of differences in the mean increment in the periods used for 
the increment index calculation (One-way ANOVA) 



Altitude 


WP00 


ANOVA 


0 

1 


ANOVA 


Ph 


p-value 


Ph 


p-value 


500 


1.01 


0.0253 


n.s. 


0.96 


0.4595 


n.s. 


600 


0.59 


4.9510 


<0.05 


1.10 


1.0891 


n.s. 


700 


0.96 


1.5399 


n.s. 


0.89 


0.2762 


n.s. 


800 


0.82 


1.4585 


n.s. 


0.80 


2.9881 


n.s. 


900 


0.81 


3.6668 


n.s. 


0.81 


5.1882 


<0.05 


1000 


1.17 


0.1744 


n.s. 


1.07 


0.5253 


n.s. 


1100 


1.14 


0.3734 


n.s. 


1.12 


0.0258 


n.s. 


1200 


1.15 


1.6167 


n.s. 


1.11 


1.1047 


n.s. 



In the phytosociological records described on the 
plots, 148 plant species were recorded: 128 species 




Fig. 4. Mean yearly increment of trees in upper forest stage at 1100 (undamaged) and 1200 (damaged by Cephalcia in 
1979-1984) m a.s.l. in the period 1961-2010 
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of vascular plants (including 26 species of trees and 
shrubs), and 20 species of mosses. Two associations and 
one community were identified (tab. 4): 

- Dentario glandulosae-Fagetum typicum Klika 1927 
em. W. Mat. 1964, represented by the typical sub-as- 
sociation in two altitudinal forms: sub-mountainous 
form and mountainous form; 

- Abieti-Piceetum Pawl, et Kulcz. 1923 em. J. Mat. 
1978; 

- Abies alba - Rubus hirtus. 



Tab. 3. Increment indices of stands in the upper forest 
zone (from 1100 m a.s.l.), corresponding to the period of 
defoliation caused by Cephalcia spp. 



Period 


1100 m a.s.l. 


1200 m a.s.l. 


•a 

Ch 


Ph 


p-value 


p-T 

$ 


PP 


p-value 


decline: 

1976-85/1966-75 


0.83 


0.5105 


n.s. 


0.71 


5.2341 


<0.05 


recovery: 

1986-95/1976-85 


0.83 


0.6868 


n.s. 


1.23 


3.9617 


n.s. 



No relationship was found between the species di- 
versity of ground vegetation, expressed by the Shannon- 
Wiener index, and level of stand diversity at individual 
altitudes (fig. 5). 




Altitude [m a.s.l.] 

Fig. 5. The Shannon-Wiener index calculated for ground 
vegetation in stands representing different tree diversity 
levels, at the various elevations 

The average Shannon-Wiener (PP) index calculated 
for pooled data concerning vegetation at four elevations 
demonstrates a clear trend of plant species diversity, 
decreasing with altitude (r 2 = 0.77, p < 0.001): the high- 
est ( H ' = 3.90) at 500 m a.s.l., the lowest ( H ' = 3.11) at 
1100 m a.s.l. (fig. 6). 



Tab. 4. Forest vegetation on 15 plots located at four altitudinal levels on Radziejowa 



Elevation 
[m a.s.l.] 


Stand 

diversity 

level 


Phytosociological plant association or community 


European forest types* 


500 


1-4 


Dentario glandulosae-Fagetum typicum sub-mountainous 
form 


7.5 Carpathian mountainous beech forest 


700 


1,2,4 


Dentario glandulosae-Fagetum typicum mountainous form 


7.5 Carpathian mountainous beech forest 


3 


Abies alba-Rubus hirtus community 


3.2 Subalpine and mountainous spruce and 
mountainous mixed spruce-silver fir -forest 


900 


1M 


Abies alba-Rubus hirtus community 


3.2 Subalpine and mountainous spruce and 
mountainous mixed spruce-silver fir -forest 


1100 


1,2 


Abieti-Piceetum 


3.2 Subalpine and mountainous spruce and 
mountainous mixed spruce-silver fir -forest 


3 


Dentario glandulosae-Fagetum typicum mountainous form 


7.5 Carpathian mountainous beech forest 



* www.eea.europa.eu. 2007. European forest types. Categories and types for sustainable forest management reporting and policy. 
EEA Technical report No 9/2006: 5—111. 
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Altitude [m a.s.l.] 

Fig. 6. Ground vegetation diversity at individual altitudes 
(level of experimental plots) 



Discussion 

The significant differences found in the value of 30- 
year radial increment of spruce populations at differ- 
ent elevations may indicate the influence of several 
factors on tree growth, other than climatic conditions. 
A decrease in the increment of spruce stands grow- 
ing in the highest parts of the Radziejowa Massif has 
already been recorded by Zawada (1994). The author 
pointed out that the defoliation caused by folivorous 
insects was the main stressor for the spruce trees that 
were weakened and then attacked by bark beetles. One 
of the measurable effects of this damage is the reduc- 
tion of radial tree growth (Capek 1962), which usually 
is temporary and - in cases where damage is moderate 
- is proceeded by rapid recovery over the following 
years. This kind of reaction, which was recorded in 
our case study but also previously shown by Jachym 
(2011) in relation to the degree of damage, suggests 
that stands on the highest parts of Radziejowa are of 
relatively high vitality and growth potential, espe- 
cially in the upper forest zone. In the Norway spruce 
stands in Western Sudetes the radial increment de- 
creased with altitude, but after an extended outbreak 
of Zeiraphera griseana (Hb.) in 1977-1983 a similar 
recovery was not recorded - probably due to a long- 
lasting water deficiency (Grodzki 1998). 

Based on the increment data, it seems that in the 
case of Radziejowa one should distinguish two main 
features of Norway spruce: productivity (expressed by 
the radial increment value), and vitality (expressed by 
the increment indices describing its dynamics). The 
variation in mean growth increment in over a 30-year 



period suggests that the optimum zone for Norway 
spruce productivity lies between 800 and 1000 m a.s.l. 
On the other hand the highest increment indices were 
found in the trees growing on higher altitudes, that in- 
dicate greater growth potential of stands in the upper 
forest stage (naturally dominated by Norway spruce), 
where the productivity of trees (with respect to incre- 
ment size) is more limited by both climatic and biotic 
factors. This was clearly demonstrated based on the 
data acquired from a relatively limited samples of trees, 
though may not be true in relation to the increment of 
stands. This result, however, concerns the growth of 
mature trees and does not apply to regeneration, as the 
optimum altitude for P. abies juveniles is lower than 
that for adults, probably because of limitations on repro- 
duction or establishment of seedlings at higher altitudes 
(Rabasa et al. 2013). 

There are many environmental factors influencing 
Norway spruce growth that differ between geographi- 
cal regions and altitudinal zones (Makinen et al. 2002). 
The lower radial increment of spruce stands growing 
on low altitude locations of the Radziejowa Massif may 
be also caused by other factors. Precipitation in the 
lower elevations is not as high as in the crest positions 
of Radziejowa (Brzezniak and Czemerda 2000), and 
fluctuating rainfall in the growing season may nega- 
tively influence the growth of Norway spruce, which 
is a species extremely susceptible to changes in water 
regime (Kurjak et al. 2010). The lack of water is the 
main growth limiting factor at low altitudes, while at 
intermediate altitudes this effect is not so strong (Maki- 
nen et al. 2002). In addition, the stands occupying lower 
elevations of the massif are characterized by a higher 
frequency of pathogenic fungi, such as white-rot fungus 
Heterobasidion parviporum (Fr.) Niemela & Korhonen, 
bulbous honey fungus Armillaria cepistipes Velen, 
and dark honey fungus Armillaria ostoyae (Romagn.) 
Herink, than those at higher altitudes of Radziejowa 
(Grodzki et al. 2012). The occurrence of these patho- 
gens is also driven by hydrothermic conditions - water 
deficiency in particular is known as a stimulating fac- 
tor for Armillaria spp. (Desprez-Lousteau et al. 2006). 
Thus, it can also be considered as one reason for lower 
growth of some of spruce stands investigated. 

The forest plant associations and communities 
found on research plots located in the Radziejowa Mas- 
sif do not represent all the variability of forest vegeta- 
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tion present. Above the range of elevations examined, 
the association Plagiothecio-Piceetum (Schaf., Pawl, 
et Kulcz. 1923) J. Mat. 1977 can be found, and below, 
the association Tilio-Carpinetum Tracz. 1962 (Mycz- 
kowski and Grabski 1962; Staszkiewicz 1995, 2000). 
Within the range of elevations studied, the association 
Luzulo luzuloidis-Fagetum Markgr. 1932. em. Meusel 
1937 also occurs, as well as communities of artificial 
origin in spruce-dominated stands, which have unde- 
termined phytosociological status (Fabijanowski 1962; 
Myczkowski and Grabski 1962; Staszkiewicz 1995, 
2000) and were widespread in the past. As studies have 
shown, despite the dominance of spruce in the stand, 
such plant communities can now be readily recognized. 
This demonstrates the significant progress of vegetation 
change processes towards sustainable forests, taking 
place in communities of planted stands where Norway 
spruce is dominant. 

The relationship between plant biodiversity and el- 
evation can be described by models (Nogues-Bravo et 
al. 2008). The main ones being: 

- a monotonic species richness model, in which spe- 
cies richness are assumed to decrease with altitude, 

- a hump-shaped species richness model, in which 
species richness is assumed to peak within the lower 
half of the gradient, 

- a mid-domain-effect model with a peak in richness 
in the middle of the gradient as a consequence of 
geometric constraints and two hard boundaries. 

The dependency of biodiversity on altitude, ob- 
served in the Radziejowa Massif, is best described by 
the first of these models. 

Conclusions 

- At the altitudes 800-1000 m a.s.l. the average radial 
increment of spruce was at its highest. It is probably 
the optimum for spruce, with favourable climatic 
conditions and low competition. The reduction in 
average increment at 900 m a.s.l. may be due to the 
occurrence of limiting factors (e.g. reduced water, 
prevalence of fungi). 

- The increase of radial increment indices with alti- 
tude and rapid recovery following insect defolia- 
tion can be considered as an indication of the high 
growth potential of spruce, which may suggest an 



altitudinal shift of the optimum growth zone for 
spruce in recent years. 

- The vegetation of Radziejowa Massif is character- 
ised by a large (148) number of species. The highest 
plant diversity, found at 500 m a.s.l., decreased with 
increasing altitude. 

- The possible increase in air temperatures may re- 
sult in changes in the vertical range of many plant 
species including trees, and as a consequence, in 
an upward shift of the boundaries of plant zones, in 
this case the sub-mountainous and lower mountain- 
ous forest zones. In this scenario, the optimal per- 
formance of Norway spruce could be limited to the 
highest altitudes. 
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